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ABSTRACT 

This work was undertaken to check the feasibility of obtaining the 
velocity profiles downstream of the rotors of turbomachines by measuring 
only the radial distributions of total temperature, total pressure, and 
absolute flow angle and using these data to calculate the flow velocities 
by the non-isentropic radial equilibrium equation such that the measured 
flow rate is matched. Curvature effects are ignored in this analysis. 

The velocities calculated by this method are compared with those 
obtained by measurements with flow probes and hot-wire anemometers. 

The tests were conducted on a 3 -stage axial flow compressor of the 
Turbo-Propulsion Laboratory at the Naval Postgraduate School. Measure- 
ments were made over the entire operating range of the compressor at a 
fixed rotor speed of about 1200 rpm. 
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SECTION 1 



INTRODUCTION 

The experimental analysis of performance of turbomachinery requires 
detailed and accurate measurement of the flow properties within the pas- 
sages of the machine under investigation. 

Even at the present state of the art it is necessary to make a num- 
ber of simplifying assumptions for the design of turbomachines. A limited 
amount of basic design data is available which may not apply directly to 
a particular type of machine . It is necessary, therefore, to check the 
accuracy of the design method by testing actual machines. For axial 
type turbomachines it is of great interest to know the flow profiles 
after the different stages since this information indicates whether the 
flow deflections, losses, and velocity distributions on which the design 
was based have been realized. Comparison of assumed and actual flow, 
properties will^make it possible to correct or improve the design pro- 
cedure with the result that the performance of later machines can be 
predicted with greater accuracy. 

The measuring of the flow properties between the rows of blades of 
axial turbomachines is connected with considerable difficulties. It is 
necessary to use small probes to avoid disturbing the flow; and some of 
the more accurate probes, such as Pitot-static probes, cannot be arranged 
because of small axial blade clearances • Flow angles, total pressures, 
and total temperatures can be measured with relative ease and sufficiently 
small probes. To establish the velocity profiles, however, it is neces- 
sary to know the static pressure and the static temperature. The meas- 
uring of these quantities between the rows of turbomachines can usually 
be carried out only with probes that have been calibrated for different 
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Mach number, Reynolds number, proximity of walls, and flow angles. 

For these reasons the accuracy of such measurements is low, and the 
experimental procedure is quite difficult. 

The purpose of this study is to determine whether it is possible 
to use the first-mentioned measurements together with certain theoreti- 
cal flow relations to obtain the flow profiles with a sufficient degree 
of accuracy. These conditions are checked on a three-stage axial com- 
pressor of large size where the axial clearances between the rows of 
blades make it possible to measure the velocities directly with cali- 
brated flow probes and hot-wire instruments. This compressor will 
henceforth be called the O.N.R. compressor since the Naval Postgraduate 
School received this machine from the Office of Naval Research, Depart- 
ment of the Navy. 

The author greatfully acknowledges the guidance and counsel of 
Professor M. H. Vavra of the Department of Aeronautics of the Naval 
Postgraduate School. Professor Vavra not only guided this research 
project but also gave the author an appreciation of what it means to be 
an engineer. 

The author also wishes to express his appreciation to the Office of 
Naval Research which made the compressor used in this study available 
to the Naval Postgraduate School. 

Thanks are due also to Mr. J. E. Hammer and Mr. R. W. Savage of 
the technical staff of the Turbo-Propulsion Laboratory as well as to 
Mr. T. B. Dunton and the staff of the Instrument Shop of the Department 
of Aeronautics. 
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SECTION 2 



BASIS FOR THE ANALYSIS 

For any flow condition there are two principal equations of conserva- 
tion which must be satisfied, namely, the equation of motion and the equa- 
tion of continuity. 

The following simplifying assumptions are made for the present treat- 
ment : 

1. The fluid is assumed to be a perfect gas. 

2. The general flow equations are applied to compute fluid proper- 
ties between blade rows where blade forces are nonexistent. 

3. The flow is assumed to be steady. 

4. The flow is assumed to be axisymmetric . By this assumption, 
all partial derivatives with respect to the peripheral coordinate are 
taken to be equal to zero. Basically, this is the equivalent of saying 
that the three dimensional flow through a blade row can be represented by 
the circumferential average of the flow conditions at each radial posi- 
tion. 

5. The flow is assumed to be adiabatic. 

6. The shearing stresses at the annulus boundaries are ignored in 
the equation of motion, but their effects are reflected by the use of 
measured data. 

7. The specific heat of the gas is assumed to be constant. 

With these simplifying assumptions, the equation of motion as defined 
by Vavra (1) and the equation of continuity, respectively, are 

VH = TVS * V * (V*V ) (1) 
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V*(<? V) = O 



( 2 ) 



Since this study deals with the radial variations of fluid properties, 
it suffices to use the radial components of Eq. 1. In consistent engi- 



neering units, 




A further simplifying assumption is made inasmuch as the radial 
accelerations due to streamline curvature are ignored, thus eliminating 
the last. term in Eq. 3. This assumption is valid for lightly loaded 
machines where the curvatures of the meridional flow path are small as 
described by Harrison (2) and Giamati and Finger (3). Since both condi- 
tions are met by the O.N*R # compressor, the error introduced by this 
assumption will be small. 

In consistent engineering units the equation of continuity can be 
expressed by 



The integration of Eq. 4 within the boundary layer of the annulus 
walls requires extensive data in the boundary layer. These data were 
not available in this study due to the physical dimensions of the probes 
which were used to obtain the surveys. The introduction of a gross 
weight flow blockage factor in the equation simulates the effect of the 
boundary layer by increasing slightly the axial velocity components to 
meet the required continuity condition. With the weight flow blockage 
factor B^, Eq. 4 becomes 



/ 




X hUb 



(4) 
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Equations 3 and 5 require iterative solutions to satisfy both equa- 
tions. In addition, the non-linearities in Eq. 3 require numerical 
solutions; hence, the use of a high speed computer is almost a necessity 
to carry out the procedure. 

Equation 3 was multiplied by dr, and integration yields 




Stations r and.r^ are two closely adjoining radial stations. 

Some of the terms in Eq„ 6 may be integrated directly, but the 

2 

integrals containing entropy changes and the one involving /r must 
be determined numerically. 

Across any blade row the difference in entropy at a given stream- 
line may be expressed by ^ 



J( s b)-J( s r) 






CE./TV .) 
(Pn/PO 



( 7 ) 



The symbol S represents the stagnation entropy. Subscript 2 refers 
to the downstream position and subscript 1 to the upstream position. 

The entropy gradient across the blade row on two closely adjoining stream- 
lines may now be calculated. For descriptive purposes two such stream- 
lines will be subscripted by i and j. 

Subtracting the axial entropy change across the blade row on stream- 
line i from that on j yields 



m,- 



iL 



f(T., t /Tr lf )*] 


TO*) 


[(ft./Pt,,) J 


L( ) - 



( 8 ) 
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The numerical value of the entropy gradient in the radial direc- 
tion may now be simplified by assuming that position 1 is far upstream 
where no radial entropy gradient exists and where radial gradients of 

total temperature and total pressure are zero. Therefore, T equals 

t li 

T and P equals P . Then Ea. 8 simplifies to 

iy hi hj 






R_i 
J ^ 






fa,/ 


* n 


k//f 


r.i) 



(9) 



The integration of the term in Eq. 6 that contains the entropy 
gradient ^ S/ c) r is then 



- i fA 



( T,/p 






Jj(T+T.) 



(10) 



The static temperature appearing in Eq. 10 may be expressed by 



the stagnation temperature T and the absolute velocity V, or by 



T = Te — 



t~l v 

2 yp 



(id 



With Eq. 11, Eq. 10 becomes 



\7?r 



I J [t /)] cm 



The non-linear term in Eq„ 6 containing /r can be evaluated 



numerically by 



n 



{ v* 



j 



ji 



d/t — 



(>'♦ ¥) 






(13) 
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For simplification, let 



F 



Jb\ 




X 



JL- 

x y-' 



(K7F7) 



(14) 



Since V at an arbitrary radius is to be determined, V and V in 
a u 

Eq„ 6 are expressed by V and o<^ , or 

a 



Vu* = Vo/ txJ 1 ©< 

V ' = Vi' ( I + ex' ) 



Equation 6 may now be solved for at an arbitrary radius r, 

= [|- F ^7- + txAo* (Z-^r-'F J 

[2,R^(T.-T« J )-F ? R(T+T*,-)-'- 

V^CF^ fH +Wx;^V.(F^!.*Z-^.)] (U) 



With measured values of total temperature, total pressure and flow 

angle as a function of radius, Eq 0 16 could be solved at all radial 

positions if the reference axial velocity ¥ were known. 

a . 

i 

In Eq* 5 the density is a function of velocity. For a perfect gas 




(17) 



From isentropic relationships 




(18) 



Combining eqs. 17 and 18 with the static temperature of Eq. 11 



This quantity introduced into Eq. 5 establishes the flow rate for 

the assumed velocity V . If it does not agree with the measured flow 

a i 

rate, a new reference axial velocity V is assumed to obtain an improved 

a . 

1 

solution of Eq. 16; and the process is repeated until agreement between 
the measured and calculated mass flow rates is reached. 

Since Eq. 16 contains differences instead of differentials, it is 
evident that small increments in radius are necessary to insure accuracy 
in the final results. 

The described method of velocity determination is similar to that 
presented in Ref. 3 for the design and off-design analysis of axial flow 
compressors . 
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SECTION 3 



COMPUTER TECHNIQUE FOR SOLUTION 

A computer program for the calculation of the radial flow properties 
downstream of a rotor row was prepared for the IBM-360-67 computer in 
Fortran IV language. 

The program is specifically set up for the 0«N.R. three-stage axial 
compressor of the Naval Postgraduate School described in Section 4, but 
slight modifications make it applicable to radial surveys downstream of 
arbitrary axial rotor rows of either turbines or compressors for the 
assumptions stated in Section 2. 

3.1 Non-Isentropic Radial Equilibrium Equation 

The computer program has been set up for data from a maximum of 
fifty radial stations, not necessarily extending to the hub and to the 
tip radii, although accuracy would be reduced if data were not measured 
close to the hub and the tip of the blading. Values of total pressure, 
total temperature, and absolute flow angle are introduced as a function 
of radius. Other input values are measured flow rate, initial reference 
value of axial velocity, assumed blockage factor, desired accuracy in flow 
rate matching, and rotor speed. 

The introduced data which are measured at discrete radii are first 
expressed with fourth order polynomials as a function of radius. The 
polynomial fit for the data was obtained by the method of least squares. 
Due to the smooth change of the data, the polynomial fit was extremely 
accurate throughout the radial span for all the parameters. An example 
of the accuracy of the curve fit is given in Table B-2 where the symbol 
X signifies radius; F2, the input data for different radii; Y, the values 
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of the parameter obtained from the polynomial; and DELY, the difference 
between the measured data and those from the polynomial. The subprogram 
used for the polynomial curve fitting is available in the Naval Post- 
graduate School subroutine library and is presented in Appendix B. 

A maximum of one hundred points may be specified at which the 
equilibrium equation is to be solved. For brevity, the non-isentropic 
radial equilibrium equation shall be referred to as the NISRE equation 
in the remainder of this presentation. The radial positions of the 
points are obtained by dividing the radius from hub to tip into equal 
increments. Values of total pressure, total temperature and absolute 
angle of flow are then obtained from the polynomials at these specified 
radial positions. The measured mass flow rate is calculated from meas- 
ured conditions in the inlet to the compressor. 

The subprogram dealing with the NISRE equation then solves for the 
axial velocity and density profile from blade hub to blade tip. The 
profiles thus obtained are used in the continuity equation which is 
integrated by the trapezoidal method. The initial axial velocity at 
the reference position is then changed arbitrarily, and a second solution 
is generated. The third and subsequent values of the reference axial 
velocity are calculated with the Newton-Raphson iteration technique. 

This process continues until the desired accuracy in mass flow rate is 
obtained. 

With the axial velocity profile thus determined, all remaining fluid 
properties may be calculated. These properties are printed out in dimen- 
sional and non-dimensional form. Examples of the print out of the solu- 
tions are presented in Tables B-3 and B-4, 
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3.2 Determination of Velocity Profiles with Flow Probes 

The measured values of total pressure, static pressure, absolute 
flow angle, total temperature and rotor speed are used to obtain the 
axial velocity profile. The fluid properties at different radial 
stations are calculated and printed out in a format similar to that of 
the NISRE equation (See Tables B-5 and B-6.). 

3.3 Determination of Velocity Profiles with the Hot-Wire Anemometer 

The hot-wire anemometer survey was made at the same positions as 

the flow probe and total temperature surveys. To place the hot-wire as 
nearly perpendicular as possible to the flow, the absolute flow angle 
at each radial position was obtained from a five-hole probe. The hot- 
wire output was set to a reference value at the mean square radius of 
the annulus, and all values were measured with respect to this setting. 
The velocity determined by the flow probe at the mean square radius was 
used as the reference velocity for the hot-wire anemometer. The root 
mean square radius was used as reference radius since the wall inter- 
ference effects on the flow probe are small at this location. 

Tables B-7 and B-8 give sample outputs of the flow properties at 
different radii calculated from the measured data. It should be noted 
that the static pressures and static temperatures are the values obtained 
from the analysis. The measured static pressures are not used in these 
calculations . 
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SECTION 4 



O.N.R. THREE-STAGE AXIAL-FLOW COMPRESSOR 

The compressor installation is shown in Fig. 1. A detailed descrip 
tion of the test installation is given by Bowen et al. (4). Marshall 
(5) describes the test set up at the Naval Postgraduate School. 

The compressor is suitable for a large variety of research experi- 
ments. Each one of the rotor and stator blade rows may be removed. 
Although stagger angles are adjustable, they were not varied during the 
present experiments. All three rotor blade cascades were set at a 
stagger angle of 43.8 degrees, with all stator blade stagger angles at 
39.8 degrees. 

The casing of the compressor was designed to give maximum accessi- 
bility for measuring instruments. Six rectangular instrument ports are 
located in the upper half of the casing. These ports accommodate a 
special instrument carriage permitting detailed surveys in radial and 
peripheral directions (Fig. 2). Three of the instrument ports are 
suitable for stator outlet surveys over a peripheral arc of 15 degrees, 
and the other three are suitable for rotor outlet surveys over the 
same arc. Only rotor outlet surveys were made during this study. 

The casing has numerous radial survey ports. Several of these 
were used downstream of each rotor. 

The compressor is electrically driven and operates at two fixed 
speeds of either about 900 rpm or 1200 rpm. 

The inlet duct (Fig. 3) consists of a screen, an entrance bellmouth 
and a length of straight pipe with a honeycomb flow straightener ahead 
of the inlet guide vanes. 
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The compressor discharge (upper foreground of Fig. 1) is fitted 
with an electrically operated throttle valve whose positions are 
indicated by a counter. 
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SECTION 5 



INSTRUMENTATION 

The compressor inlet duct allows vertical and horizontal surveys 
to be made across the inlet as well as at 60 degrees from the vertical 
direction by unbolting and rotating the inlet duct. A 3/32 inch diam- 
eter Prandtl-type Pitot-static probe was mounted for vertical survey 
(Fig. 4) with readout on a Meriam Instrument Company Model-750 Micro- 
manometer (Fig. 5) • A 1/4 inch diameter hot-wire anemometer probe was 
mounted for horizontal survey. The hot-wire probe was constructed at 
the Naval Postgraduate School and provides for quick changing of tips 
in the event of wire breakage (Fig. 6) „ A supply of tips was kept 
readily available during the inlet surveys, and tip replacement was 
effected on several occasions while the compressor was in operation. 
Hot-wire data were read on a Security Associates Applied Sciences Model 
200 Constant Temperature Anemometer. Details of the flow rate measure- 
ment are given in Section 6. 

Wet and dry bulb Mercury glass thermometers were mounted on the 
inlet screen. 

A 1/4 inch diameter Prandtl-type Pitot-static probe was installed 
between two of the inlet supports ahead of the guide vanes (foreground 
of Fig. 7) with readout on the same micromanometer as the inlet duct 
Pitot-static probe. This was made possible by the use of two scanner 
valves (Fig. 8). In addition, a 1/8 inch diameter Kiel probe with a 
total type thermocouple was mounted in the same meridional plane through 
one of the radial access ports. This total temperature probe is shown 
in the foreground of Fig. 7. The reference junction of the thermocouple 
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was placed in an ice bath. The output of the thermocouple was read on 
a Leeds and Northrup Potentiometer (Fig. 9) . 

Radial surveys were conducted downstream of each of the three rotor 
blade rows. Total and static pressure measurements were obtained with 
a United Sensor Corporation DA120 five hole probe (serial number 535) 
although the pitch angle was assumed to be zero. The five hole probe 
was mounted in the special carriage mentioned in Section 4 with readout 
on a water manometer board (Fig. 10). Flow direction was obtained by 
balancing pressure on the static ports. Total pressure and the differ- 
ence between total pressure and static pressure were measured with the 
same manometer board. Total temperature measurements at different radii 
were obtained by use of a thermocouple mounted in a 1/4 inch diameter 
Kiel probe which was placed in one of the radial access holes. The 
reference junction for this thermocouple was also placed in an ice bath, 
and the aforementioned potentiometer was used for the readout by the use 
of the double-pole double-throw knife switch shown in the background of 
Fig. 9. Hot-wire anemometer data at different radii were measured after 
each of the three rotor blade rows by means of a 1/4 inch diameter hot- 
wire probe. Figure 11 shows the installation of the total temperature 
probe and the hot-wire anemometer probe used for the radial surveys. The 
hot-wire anemometer readout instrument is shown in Fig. 9. 

Compressor speed was measured with a flux cutter mounted on the 
compressor drive shaft whose impulses were read on a Monsanto Model 100A 
Counter-timer. 
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SECTION 6 



FLOW RATE CALIBRATION 

The purpose of the flow calibrations was to obtain a factor which, 
when multiplied by the velocity ahead of the inlet guide vanes, gives 
the volume flow rate through the compressor., 

Marshall (5) found a flow rate calibration factor for the O.N«R. 
compressor by means of horizontal and vertical surveys of the inlet with 
the existing Pitot-static probe. This calibration factor was to be 
checked by a more accurate survey with a hot-wire anemometer. 

Since Marshall showed that the velocity profiles in the vertical 
and horizontal planes were similar, it was decided to conduct a survey 
in the horizontal plane only. 

Since hot-wire instruments give only the ratio of the actual velocity 
to some reference velocity, initial hot-wire instrument setting was made 
at the inlet duct centerline where the velocity obtained from readings of 
the Prandtl type Pitot-static probe was used as the reference velocity. 

At the inlet duct centerline the Pitot-static probe is 18 inches from the 
duct wall; hence, wall effects are minimized * The hot-wire instrument 
was set at 90 per cent of full readout range at the centerline, and 
surveys were made at radial increments of 3 inches to a radius of 15 
inches and of 0.2 inches to a point 0*2 inches from the duct wall. The 
design of the hot-wire probe prevented measurementcS closer to the wall. 
Frequent checks were made at the duct centerline to ascertain that the 
centerline setting of 90 per cent had not drifted. Such instrument 
drifts were noticed only on rare occasions. If drift occurred, it was 
always less than 2 per cent of full scale, and the hot-wire instrument 
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was then reset to 90 per cent. Surveys of this type were conducted on 
both sides of the centerline. Total pressure and velocity head readings 
were also taken at the duct centerline and ahead of the inlet guide 
vanes with the fixed Pitot-static pressure probe connected to a water 
micromanometer. Wet and dry bulb temperatures were recorded from the 
Mercury thermometers on the inlet screen. The temperature in the test 
cell was recorded also for the purpose of correcting water column 
heights . 

The following symbols are used for the recorded data: 



q lr 


-velocity head at the centerline of the inlet duct (in.H^O) 


P tl 

lr 

q 2r 


-total pressure at the centerline of the inlet duct (in.H 9 0) 
-velocity head ahead of the inlet guide vanes (in. 11^0 ) 


TOR 


-total pressure ahead of the inlet guide vanes (in. H^O) 
-temperature at the inlet screen (°F) 


BAR 


-barometric pressure (in.Hg.) 


TBAR 


-temperature of the Mercury column (°F) 


TW 


-temperature of the water column (®F) 


VRR 


-hot-wire anemometer readings, right of centerline looking 

downstream (dimensionless) 


VRL 


-hot-wire anemometer reading, left of centerline looking 
downstream (dimensionless) 



The specific gravity of the Mercury in the barometer was corrected 
for temperature by 

GHGRM = l3.U305’-0.00l3U03O3TBA*R <20) 
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The specific gravity of the water column was corrected for tempera- 
ture by 

GWRn = O.S9837£33+|.0605V5£TW‘IO' 4 -!.S'93f86Tw‘* 10“ 6 < 21 > 

These specific gravity relations were obtained from tabulated data 
of Ref. 6. 

With these corrections, pressures were converted to absolute values, 
or 

PAMB = BAR C0.4S3lf85")(;H6RM ( psfe.) (22) 

where the constant 0.4891585 is the conversion factor from inches of 

Hg. to pounds per square inch for a specific gravity of Mercury of 13.54. 

Then 

a,= Pt.J&WRM) & z -^ 73g +PAI-1B (psfO < 23 ) 

The constant 62.42732 is the conversion factor from feet of water 
to pounds per square foot for a specific gravity of water of 1.0. 

%. = t ,«( GvRM) a '^ 792 (psf) <24) 

The static pressure was then calculated by 

P| = Pt- ( psf A,) (25) 

Because of the low velocities of about 70 feet per second in the 
inlet duct, and ahead of the compressor inlet guide vanes, compressi- 
bility effects were ignored, and the temperature TOR at the inlet screen 
was used as the static temperature * 

T 0 = TOR + 459.7 CR) < 26 > 
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Considering the air entering the compressor as a perfect gas having 



a gas constant R of 53.35 (ft-lb)/(lbm-°R) , the mass density is 



p _ -£l_ 



(slug/fi 3 ) 



(27) 



For each run the density was taken as the average obtained from all 
readings . 

The velocity at the inlet duct centerline was then calculated by 



V, 



2ft, 

1 °. 



(f t/s&c ) 



(28) 



This velocity was used as the reference velocity for the hot-wire 
anemometer setting. Then, the average velocity V at an arbitrary radius 
is taken as 



V = (m i . 2 VRR) ^ (ft/4et) 



(29) 



The constant 0.9 must be introduced since the hot-wire anemometer 
was set at 90 per cent of full scale readout of the instrument at the 
centerline in case velocities higher than those at the centerline would 
occur. 

The volume flow rate was calculated by 



V0l_FL= 3.1416 V, (4rf+ 3.1416 zivj ( * ■ *,+**- ; ) * ~ 

2 * 






;t3.l4l6 V.fcO*- (ft 'Me) «« 



The subscript 1 in the summation refers to the inlet duct center- 
line. The subscript n pertains to the reading 0.2 inches from the duct 
wall, and the subscript n+1 is for the conditions at the duct wall. 
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The velocity V n ahead of the inlet guide vanes was obtained in a 
similar manner, using the temperature measured at the inlet screen as 
the static temperature since the difference between total and static 
temperature is less than 0.4®F at the highest velocity of 70 feet per 
second. 

The calibration constant for each run was obtained from 

CCC =• — -9 !r - (ftVscc per ft/sec) ( 31 ) 

3 

This constant was found to be 4.4050 ft /sec per ft/sec velocity 
obtained ahead of the inlet guide vanes. The maximum relative devia- 
tion of CCC for all surveys was 0.3 per cent. 

This data reduction was carried out on the IBM-360-67 computer, 
program is explained in Appendix A which also gives a sample of the 
computer output. 



The 
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SECTION 7 



RADIAL SURVEYS 

Radial surveys were conducted downstream of each row of rotor blades 
of the O.N.R, compressor. Radial distributions of total pressure, veloc- 
ity head, total temperature, absolute flow angle, and hot-wire velocity 
ratio were obtained at various flow rates. 

Flow rate was varied by means of the discharge throttle valve. After 
setting the throttle valve to a new position, data were taken only after 
waiting for about 20 minutes to stabilize the temperature in the test 
cell. The test cell has a front door through which air is normally admit- 
ted to the compressor from the atmosphere. The air is discharged back 
into the atmosphere through a door located opposite to the front door and 
aft of the compressor. However, experience showed that gusts in the 
atmosphere produced erratic measurements particularly with the hot-wire 
instrumentation. By closing the front door it was found that the flow 
oscillations could be minimized, but the above-mentioned waiting period 
was necessary for the air to reach a new equilibrium temperature in the 
test cell. Steady temperature conditions were checked by the mercury 
thermometer at the inlet screen and the total temperature Kiel probe 
ahead of the inlet guide vanes. 

Periodic readings of the total pressure, velocity head and total 
temperature ahead of the inlet guide vanes were recorded during each run 
for the determination of the flow rates. Conditions at the inlet screen 
and ahead of the inlet guide vanes were found to be very stable for all 
flow rates irrespective of throttle settings and atmospheric conditions. 

The previously described traversing instruments for the radial 



35 



surveys were mounted downstream of each one of the rotors. Initially, 
a peripheral survey was made with the flow probe at various radii to 
check whether the flow conditions were uniform in the peripheral direc- 
tion. Generally, this was found to be true downstream of each row of 
rotor blades at flow rates above full stall. 

For most of the radial survey runs a hot-wire instrument setting of 
90 per cent was used at a reference radius of 14 . 5 inches, or very close 
to the root mean square radius of the annulus. To properly position the 
hot-wire with respect to flow direction, initial angle setting of the 
probe was made with the flow angles obtained from the flow probe. This 
setting was then refined by turning the hot-wire into the direction 
which gave maximum output. Almost without exception, the hot-wire angle 
position thus found was identical with the angle measured by the flow 
probe. The instrument setting of 90 per cent at the radius of 14.5 
inches was periodically checked for drift. Drift was encountered only 
occasionally and was always less than 2 per cent. If drift was noticed, 
the hot-wire instrument was reset to 90 per cent of full scale at the 
reference radius. 

The total temperature Kiel probe was set at the angle which was 
obtained from the flow probe at the particular radius. 

For each run, radial survey readings of total temperature, total 
pressure, velocity head, absolute flow angle, and hot-wire anemometer 
velocity ratio were taken from the blade tip radius of 18 inches to a 
distance of 0.2 inches from the hub which has a radius of 10.8 inches. 
Although readings closer to the hub would have been desirable, the design 
of the flow probe prevented readings at distances less than 0.2 inches 
from the hub. 
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Data were taken at radial increments of 0.5 inches for the majority 
of the runs. Initially, the increment in radius was chosen as 0.2 inches. 
However, because of the smooth variation of the flow conditions in the 
annulus, this increment could be increased to 0.5 inches without loss 
of accuracy. 

Various flow rates were set by means of the discharge throttle valve. 
The throttle valve position is indicated by means of a digital counter. 

The counter is calibrated from 0 to 1360 (closed to open). Radial surveys 
downstream of each rotor row were taken at counter readings of 1000, 900, 
800, 700, 600, 500, and 435 - Below a counter reading of 435, low rumbling 
noises and vibrations were observed as well as significant read-out fluc- 
tuations, indicating that major portions of the blading were stalled. 
Peripheral surveys with the flow probe at a counter setting below 435 
showed significant flow variations with respect to the peripheral direc- 
tion. These conditions existed between the hub and approximately mid- 
blade height. In addition, extreme variation in the hot-wire anemometer 
reading was noticed also in this region. In some cases, the hot-wire 
oscillations approached 30 per cent of the instrument read-out. Because 
of this condition no surveys were made at throttle settings below 435 
counts. At a throttle setting of 435 counts, slight fluctuations in the 
hot-wire anemometer readings were noticed between the hub radius of 10.8 
inches and a radius of 12 inches, whereas the flow probe readings were 
steady in this portion of the annulus. 

No flow profiles are presented for readings of the throttle valve 
counter above 1000, although four such surveys were undertaken. They do 
not show appreciable differences since the flow rate increases only by 
about 2 per cent for counter readings between 1000 and 1360. 
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SECTION 8 



DISCUSSION OF RESULTS 



8.1 General 

For the present study the 0*N.R. compressor was run for 68 hours, 
of which 28 hours were required for the flow rate calibrations . 

The tests were conducted with the so-called free-vortex blading 
of the compressor installed. The rotor blades of all three rotor rows 
were set at the design stagger angles of 43.8 degrees with respect to 
the axial direction. All stator blades had stagger angles of 39.8 
degrees. The runs were made at the high rotor speed of about 1190 rpm. 

8.2 Flow Rate Calibration 

The flow rate calibrations were undertaken to verify the calibra- 
tion factor found by Marshall (5), The fluctuating readings due to 
atmospheric conditions experienced by Marshall were eliminated by clos- 
ing the door in front of the compressor inlet duct. It is interesting 
to note that during the radial surveys downstream of the third rotor, 
the door to the atmosphere was opened. Even though the atmospheric 
conditions were calm, fluctuations of up to 10 per cent were noticed 
in the hot-wire readings. The fluctuations stopped when the door was 
closed. 

Velocity measurements in the inlet duct were made from the center- 
line to a radius of 17.8 inches. It was not possible to measure the 
velocities closer than 0.2 inches from the wall. Since the inlet pipe 
has a diameter of 36 inches, the hot-wire probe had to be supported 
from both walls to eliminate vibration. Therefore, the probe had to 
be attached at right angles to the support structure (Fig. 6). The 
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probe thickness of 1/4 inch and the fillet at the attachment point 
prevented measurements closer than 0.2 inches from the wall. Velocity 
measurements to the wall would have been desirable for a more accurate 
calculation of volume flow rate. In the calculation of volume flow 
rate, the velocity measured at a radius of 17.8 inches was assumed to be 
that at the wall also. 



8.3 Presentation of Data 

The axial velocity profiles obtained from the NISRE solution and 
from the flow probe and hot-wire measurements are presented in Figs. 

12 through 32 in dimensionless form. These quantities were made dimen- 
sionless by dividing the velocities by the rotor tip speed. The parame- 
ter (p thus obtained is convenient for two reasons. First, assuming 

that a mean axial velocity V exists in the flow annulus, the volume 

a 

f 1 ow rate can be expressed by 

vr= \£ -^rar- (-ftVsec) (321 



and the rotor tip speed is 

U T - Cx> 



( f t / see) 



(33) 



where CO is the angular velocity of the rotor in radians per second. 
Thus 



m = V£ 

V CO 



ii y )44- \ _ 

18A65-I.4-/ 



CO 



C 0. 1 4-74-) 



(34) 



Therefore, (p is a flow coefficient relating volume flow rate to 
rotor speed. A numerical average of the values of CP obtained from 
the radial distribution of axial velocities at a given flow rate gives 
an average value of (P which is a measure of the volume flow rate at 
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a particular throttle setting and rotor speed. Second, the flow coef- 
ficient is a convenient parameter for comparing the velocities obtained 
by the NISRE solution and those measured with the flow probe and the 
hot-wire anemometer because changes in rotor speed have no influence on 

<p. 

The total pressures measured with the flow probe downstream of 
each of the rotors at various flow rates are presented in Figs 33 through 
35. These pressures are shown in dimensionless form by the parameter 

V = <35 > 

o 

where the subscript "a® 1 refers to the atmospheric conditions. 

The work input of a blading element in an axial-flow compressor 
can be expressed by the so-called Euler turbine equation. The rise in 
total temperature from station (1) ahead of to station (2) after the 
blading is 

&T ( = c^^Vu- VO CR) <36, 

For a frictionless process of an incompressible fluid, the theoretical 
rise in total pressure ( A then becomes 

(APe)*= pcoOtjV^-JT.V'u.,) (ps-f) (37) 

With losses, the pressure rise /S is 

APr = n = ( OW ? 0t,Vu l -K.V«.) (pi-f) (38) 

where ^ is the stage efficiency. Rewritten 

AP t = PUt’I = pU T ‘^ 
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is therefore a function of the 



The so-called pressure coefficient 
flow angles and the losses only and is independent of the density of 
the flow and the rotor speed, as long as the influences of Reynolds number 
and Mach number can be ignored. The aforementioned parameter Ip is suit- 
able for comparing the radial distributions of the total and static pres- 
sures measured after the rotors because the influences of varying densi- 
ties and rotor speeds are eliminated. 

Since the present study was not concerned with determining blade 
losses, no measurements of total pressure were made upstream of the rotors 
while the radial surveys were conducted downstream of them. Since the 
surveys after the rotors were conducted on different days, it was not 
possible to maintain exactly the same value of the average flow coef- 
ficient CP ; hence, the values of 'Ip presented for the total pressure 






cannot be used to establish the quantity vy defined earlier. 

The static pressures calculated by the NISRE method, those measured 
by the flow probe, and the ones calculated by using the velocity measured 
by tl)e hot-wire anemometer were made dimensionless by the same method. 

The dimensionless parameter thus obtained is convenient for comparing the 
static pressures determined by these different methods because changes of 
atmospheric pressure, density, and tip speed do not influence the data. 

Figures 36 through 56 show the measured values of the absolute flow 
angles 0( obtained for all runs. 



8.4 Accuracy of the NISRE Method 

Figures 12 through 32 show very good agreement between the velocity 
profiles calculated by the NISRE method and those measured by the flow 
probe and hot-wire anemometer. The agreement is good downstream of all 
rotors through the range of flow rates, with the possible exception of 
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Fig. 24. This figure shows axial velocity profiles obtained after the 
second stage rotor at a flow rate of 22.85 lbm/sec corresponding to a 
discharge throttle valve counter setting of 500. This run was repeated 
to check the results: but for both runs they were almost identical. 

The calculated velocities differ from the measured ones by as much as 

7.5 per cent at radius ratios R/R^ of 0.67 and 0.9. No obvious explana- 
tion could be found for this deviation from the otherwise close agree- 
ment between the calculated and measured data. 

Figures 57 through 77 show the comparison of the static pressures 
obtained from the three methods. This is not an independent comparison 
since these pressures are dependent on the velocities also. With the 
flow probe, the static pressures were measured directly, but for the 
NISRE calculation and hot-wire measurement, the static pressures were 
calculated by using Eqs. 11, 17, and 19. However, good agreement has 
been obtained with the three methods. 

8.5 Blockage Factor 

The value of the gross weight flow blockage factor to be used with 
the NISRE method has a large influence on the results. Reference 2 
suggests a blockage factor of 0.98 for the first stage and 0.96 for all 
subsequent stages. 

Blockage factors of 0.96 for the second and third stage rotors gave 
the results of Figs. 19 through 32. Blockage factors of 0.98 were first 
tried for the calculations of the velocity profiles for the first stage 
rotor but gave velocities well below those measured by the flow probe and 
the hot-wire anemometer. The effect of the area blockage by the wall 
boundary layers was then increased by changing the blockage factor to 0.97 
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and finally to 0.96. This increase in area blockage did not materially 
change the shape of the calculated profiles but increased the magnitudes 
of the velocities. Blockage factors of 0.96 for the first rotor gave the 
best correlation between the axial velocities calculated by the NISRE 
method and those measured. Figures 12 through 18 show the axial velocity 
profiles after the first rotor for this value. It is believed that the 
requirement for a smaller blockage factor after the first rotor is caused 
by the boundary layer build-up along the walls of the inlet duct which 
has a length of six feet. 

8.6 Compressor Design Point 

As stated in subsection 8.1, the tests were conducted with the so- 
called free-vortex blading of the compressor. At their design point, such 
bladings are supposed to produce flow patterns where the product of the 
radius r and the peripheral component of the absolute velocity remain 
constant in the radial direction after each row of blades. If entropy 
gradients are ignored, free-vortex flows have constant total pressure P , 
constant total temperature T ^ , and constant axial velocity components 
between the rows of blades. However, if entropy gradients occur even 
though total temperatures do not change in the radial direction, the axial 
velocity components of the absolute velocity are not constant. Equation 
3, without the term involving meridional streamline curvature, can be 
expressed also by 
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Hence, if the product rV^ were constant in the radial direction, the 
term containing would be zero. If the total temperature were also 
constant in the radial direction, the only terms remaining are those 
involving entropy and axial velocity component. This indicates that 
with the aforementioned assumptions, the variations of the axial velocity 
components are a function of the entropy gradients and the static tempera 
tures. 

From the presented data it can be noticed that neither constant 

velocities V nor constant total pressures P were reached with the 
a t 

installed blading, irrespective of flow rate. In Ref. 3, the design 

stagger angle of the stator blades is given as 29.8 degrees. Hence, 

since a value of 39.8 degrees was used for the present tests, it is 

likely that radial changes in V and P were due to this condition. 

a t 

All the axial velocity profiles calculated by the NISRE method show 
regions of increased axial velocities at radius ratios R/R^ °f about 
0.67 and 0.90. At which of these two stations the maximum velocity 
occurs depends on the flow rate. As flow rate decreases, the region of 
highest velocity moves from the inner to the outer part of the annulus. 
The measured velocity profiles show the same trend but are generally 
flatter. The crossover point, where the maximum velocity at both regions 
has about the same magnitude, occurs at a flow rate of about 26.5 lbm/sec 
as shown in Figs. 12, 14, 15, 22, and 29. To this flow rate corresponds 
a flow coefficient Cp of about 0.44, which is the average of the flow 
coefficients obtained from the calculated and measured axial velocity 
profiles presented in the abovementioned figures. Except for the regions 
near the hub and the tip, the axial velocity is nearly constant at this 
flow rate. Further evidence that the design criteria of a free-vortex 
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flow were most nearly met at a (1/ of 0,44 is given in Figs. 33 through 
35, where it is shown that total pressure is nearly constant except for 
the hub and tip regions. However, this uniformity is not so well real- 
ized after the second and third stage rotor. For this discussion it is 
of interest to note that Ref. 4 gives a design value of of 0.45 for 

the free-vortex blading at rotor and stator stagger angles of 43.8 and 
29.8 degrees respectively. 

Reduced total pressures and low axial velocities occur in the regions 
near the hub and near the outer casing downstream of all rotor rows. 

They are caused by end and tip clearance losses, by the wall boundary 
layers, and by the effect of the rotating hub. These conditions produce 
large entropy gradients at the inner and outer radii of the annulus for 
the reason shown in the explanation of Eq. 3 in this subsection. Their 
effect near the hub extends from radius ratios of 0.6 at the hub to 
about 0.64 for the first rotor, 0.66 for the second rotor, and 0.69 for 
the third rotor. Near the outer casing, they influence the velocity 
distributions at radius ratios of about 0.91 after the first rotor and 
about 0*94 after the second and third rotor. The regions which are af- 
fected by the hub and tip losses become larger as the fluid proceeds 
downstream since the losses accumulate from stage to stage and further 
deteriorate the flow profiles. 

The velocities at the outer casing are difficult to measure. Since 
the probes extended into the flow annulus from the outer casing, radial 
position changes of the probes of the order of 0.01 inches at the outer 
radius caused changes in the measured data of up to 40 per cent. It is 
believed that these discrepancies were caused by small errors in the probe 
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position indication. The probes retract into the radial access holes, 
and small changes in radial position cause the access holes to block 
the measuring holes of the probes. 
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SECTION 9 



CONCLUSIONS AND RECOMMENDATIONS 

The accurate determination of the radial distribution of absolute 
velocity in magnitude and direction is important for the evaluation of 
the performance of turbomachines. If these distributions were known, the 
velocity triangles would be completely determined. Knowledge of the veloc- 
ity triangles allows direct comparison between the applied design criteria 
and the actual results for a turbomachine. 

The accurate predictions of these velocity triangles is of particular 
importance in the design of multi-stage axial-flow machines. With the 
exception of the first stage, the inlet conditions of a given stage are 
the discharge conditions of the preceding stage. If a stage has inlet 
conditions different from those for which it was designed, its perform- 
ance is likely to be marginal and could produce discharge conditions 
that deviate even more from the design values. The result is that the 
machine will not be able to handle the desired enthalpy change at the 
expected level of efficiency. The radial distributions of losses and 
flow deflections in a row of blades are a function of inlet conditions 
as well as the design of the blades. A blade element designed for a 
certain range of flow incidence angles, but operating outside of this 
range, may either stall or not develop the required blade forces, de- 
pending on whether the incidence angle is too large or too small. Optimum 
design of multi-stage axial-flow machines therefore requires that the 
discharge conditions of any stage match the design inlet conditions of 
the following stage. 

Since designs of turbomachines are based on the results of experimental 
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work and the experience of the designer, checking the results of new 
designs is important for two reasons: first, to evaluate the perform- 

ance of a new design, and second, to improve the design procedures for 
future machines. 

The NISRE method represents a fairly accurate and simple method for 
determining the velocity profile downstream of an axial-flow rotor. The 
data which are required can be obtained by small and accurate probes in 
the vaneless space between blade rows with small axial clearance. The 
method as presented in this study is limited by the fact that the ef- 
fect of meridional streamline curvature is ignored. For designs where 
the curvatures of the meridional streamlines are appreciable , the comput- 
ing process could be improved in the manner carried out by Harrison (2). 

Further experimentation seems necessary to accurately predict the 
gross weight flow blockage factors, particularly with machines having 
many stages. The blockage factors used for the velocity profiles after 
the first rotor were determined by arbitrary variation until accurate 
matching with the measured profiles were obtained. By changing the 
length of the inlet duct of the O.N.R* compressor the effect of boundary 
layer growth on blockage factors could be investigated. 

The hot-wire anemometer was not used for an analysis of the unsteady 
flow in the passages downstream of the rotors. The presented average 
velocities were read on an instrument which heavily damps the signal and 
cannot follow the high frequency response of the hot-wire anemometer 
which is up to 50 kilohertz. However, some oscilloscope traces were 
photographed downstream of the second rotor. Sample pictures of these 
oscillograms, with the hot-wire probe at a radius of 13 inches, are pres- 
ented in Fig. 78. The oscillogram marked (a) was made with the throttle 
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valve open, and the oscillogram marked (g) was made at a throttle valve 
counter setting of 425. The measured flow rates are indicated in the 
figures. The time for successive blades to pass by the hot-wire is 
1.675 milliseconds, corresponding to a frequency of about 0.6 kilohertz 
which is well within the frequency response of the hot-wire anemometer. 
Therefore, the hot-wire probe will be able to measure the flow condi- 
tions in the wakes between neighboring blades. 

The low signal amplitudes in Fig. 78, which are marked by arrows, 
were found to have the time interval of 1.675 milliseconds. Hence, they 
are caused by blade wakes. The trace of the hot-wire output moves from 
left to right on the oscillograms. Hence, the wakes start on the con- 
cave side of the blade profiles and end on the convex side of the pro- 
files. Therefore, the flow profiles from left to right between the 
indicated wakes start from the pressure side of one blade and extend to 
the suction side of the neighboring blade. 

An interesting phenomenon occurs in the oscilloscope traces at a 
flow rate of 19.26 lbm/sec which is evident from the oscillogram marked 
(g). An additional trace of lesser intensity but greater amplitude is 
superimposed on the aforementioned wake pattern. This pattern is periodic 
in nature and seems to be an indication of the propagating stall discussed 
by Iura and Rannie (7). 

If the resolution of the oscilloscope traces could be improved or if 
the hot-wire signal were recorded on magnetic tape and transferred to a 
strip chart recorder, the wake pattern of the flow profiles could be 
analyzed in greater detail to detect incipient stall conditions. Using 
instrumentation of this type, it seems feasible to apply the momentum 
theorem to determine the losses associated with either different blade 
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shapes or various stagger angles of a given set of blades. The O.N.R. 
compressor is suitable for such studies since it is possible to vary 
rotor stagger angles and to install different types of blades. Such a 
study would require accurate calibrations of the hot-wire probes to 
enable measurements of velocity in absolute values rather than as ratios 
with a known velocity. A test rig is now being designed at the Naval 
Postgraduate School for such calibrations. 
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FIGURE 3 

COMPRESSOR INLET DUCT 
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FIGURE 4 

INLET WITH INSTALLED FLOW RATE SURVEY INSTRUMENTS 
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FIGURE 5 
MICROMANOMETER 
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FIGURE 6 

DETAIL OF INLET HOT-WIRE ANEMOMETER PROBE 
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FIGURE 9 

HOT-WIRE ANEMOMETER INSTRUMENT READOUT 
AND THERMOCOUPLE READOUT 
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VI MW OF COMPRESSOR SHOW TNG FLOW PROBE CARRIAGE 
INSTALLED AND ASSOCIATED MANOMETER BOARD 
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FIGURE 11 

VIEW OF COMPRESSOR SHOWING HOT-WIRE ANEMOMETER 
AND TOTAL TEMPERATURE PROBES INSTALLED 
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FIGURE 16 

AXIAL VELOCITY RATIO VS. 
RADIUS RATIO 
FIRST STAGE ROTOR 
W = 24.52 lbm/sec 
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FIGURE 17 

AXIAL VELOCITY RATIO VS. 
RADIUS RATIO 
FIRST STAGE ROTOR 
W = 22.35 lbm/sec 
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FIGURE 34 

TOTAL PRESSURE COEFFICIENT VS. RADIUS RATIO 
SECOND STAGE ROTOR 
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FIGURE 35 

TOTAL PRESSURE COEFFICIENT VS. RADIUS RATIO 
THIRD STAGE ROTOR 
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FIGURE 36 

<X VS. RADIUS RATIO 
FIRST STAGE ROTOR 
W = 28.49 Ibm/sec 




FIGURE 38 

ex' VS. RADIUS RATIO 
FIRST STAGE ROTOR 
W = 26.97 lbm/sec 
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FIGURE 40 

ex' VS. RADIUS RATIO 
FIRST STAGE ROTOR 
W = 24.52 lbm/sec 




FIGURE 42 

O' VS. RADIUS RATIO 
FIRST STAGE ROTOR 
W = 20.75 lbm/sec 




FIGURE 41 

c< VS. RADIUS RATIO 
FIRST STAGE ROTOR 
W = 22.35 lbm/sec 
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FIGURE 43 

o< VS. RADIUS RATIO 
SECOND STAGE ROTOR 
W = 28.79 lbm/sec 



FIGURE 44 

c* VS. RADIUS RATIO 
SECOND STAGE ROTOR 
W = 28.09 lbm/sec 





FIGURE 45 

ex’ VS. RADIUS RATIO 
SECOND STAGE ROTOR 
W = 27.16 lbm/sec 



FIGURE 46 

o< VS. RADIUS RATIO 
SECOND STAGE ROTOR 
W = 26.46 lbm/sec 
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CX VS. RADIUS RATIO 
SECOND STAGE ROTOR 
W = 24.83 lbm/sec 




FIGURE 48 

o< VS. RADIUS RATIO 
SECOND STAGE ROTOR 
W = 22.85 lbm/sec 




FIGURE 49 

cX VS. RADIUS RATIO 
SECOND STAGE ROTOR 
W = 20.36 lbm/sec 
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o< VS. RADIUS RATIO 
THIRD STAGE ROTOR 
W = 28.62 lbm/sec 




FIGURE 51 

o< VS. RADIUS RATIO 
THIRD STAGE ROTOR 
W = 27.91 lbm/sec 




FIGURE 52 

c< VS. RADIUS RATIO 
THIRD STAGE ROTOR 
W = 27.19 lbm/sec 




FIGURE 53 

ex' VS. RADIUS RATIO 
THIRD STAGE ROTOR 
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FIGURE 55 

o' VS. RADIUS RATIO 
THIRD STAGE ROTOR 
W = 22.05 lbm/sec 
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FIGURE 56 

oC VS. RADIUS RATIO 
THIRD STAGE ROTOR 
W = 20.71 lbm/sec 
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FIGURE 57 

STATIC PRESSURE COEFFICIENT 
VS. RADIUS RATIO 
FIRST STAGE ROTOR 
W = 28.49 lbm/sec 
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FIGURE 60 

STATIC PRESSURE COEFFICIENT 
VS. RADIUS RATIO 
FIRST STAGE ROTOR 
• W = 26.17 lbm/sec 



<N 



0 



Q. 

< 

a. 

i 

ft. 



- 0.1 



NISRE 

O FLOW 
□ HOT-WIRE 




- 0.2 




0.7 R/R Q 0.8 



0.9 



1.0 



84 



□ O 



P-P A /yOU, 





85 



n<y/ v a-a 




86 



0.2 



FIGURE 64 

STATIC PRESSURE COEFFICIENT 
VS. RADIUS RATIO 
SECOND STAGE ROTOR 
W = 28.79 lbm/sec 
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FIGURE 66 

STATIC PRESSURE COEFFICIENT 
VS. RADIUS RATIO 
SECOND STAGE ROTOR 
W = 27.16 lbm/sec 
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FIGURE 68 

STATIC PRESSURE COEFFICIENT 
VS. RADIUS RATIO 
SECOND STAGE ROTOR 
W = 24.83 lbm/sec 
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FIGURE 77 

STATIC PRESSURE COEFFICIENT 
VS. RADIUS RATIO 
THIRD STAGE ROTOR 
W = 20.71 lbm/sec 



0 . 



1 



0.6 



0.7 R/R q 0.8 



0.9 1.0 



94 





1 . 67 5 cm 



(b) W 



2 7.91 



cm 



Lbm/ sec 



Ltuili 1 1 1| 

0 1 2 



PROBE POSITION OSCILLOSCOPE SETTINGS 

DOWNSTREAM OF TIME SCALE 1 millisec/cm 

SECOND STAGE ROTOR AMPLITUDE 10 millivolt/cm 

RADIUS =- 13 inches 

HOT-WIRE CURRENT 25 milliamps 

FIGURE 7b 

HOT-WIRE ANEMOMETER OSCILLOGRAMS 
SECOND STAGE ROTOR 

‘>5 







1.675 cm 



(c) W = 27.12 Ibm/sec 





96 










1.675 cm — H 



(f) W = 22.52 lbm/sec 




( g ) W = 19.26 lbm/sec 0 1 2 



FIGURE 78 (continued) 



‘17 





REFERENCES 



1. Vavra, M. H. , Aero-Thermodynamics and Flow in Turbomachines , 

John Wiley & Sons, Inc., New York, N. Y.,1960. 

2. Harrison, Robert Glen, "An Analysis of Single Stage Axial-Flow 
Turbine Performance Using Three-Dimensional Calculating Methods," 
Naval Postgraduate School Thesis, September 1967. 

3. Giamati, Charles C. Jr. and Finger, Harold B. , "Design Velocity 
Distribution in Meridional Plane," Chapter VIII of NASA SP-36 , 
Washington, D. C. , 1965. 

4. Bowen, John T., Sabersky, Rolf H. and Rannie, Duncan W., 
"Theoretical and Experimental Investigations of Axial Flow 
Compressors (Summary Report)," California Institute of Technology, 
January 1949. 

5. Marshall, Bruce Cameron, "Effect of Stator Blade Orientation on 
the Performance of an Axial Flow Compressor," Naval Postgraduate 
School Thesis, September 1967. 

6. National Research Council, International Critical Tables , McGraw- 
Hill, 1928. 

7. lura, T. and Rannie, W. D., "Observations of Propagating Stall in 
Axial-Flow Compressors," California Institute of Technology, 

April 1953. 



98 



APPENDIX A 



DETAILS OF FLOW RATE CALIBRATION COMPUTER PROGRAM 

This program uses Che data obtained from the hot-wire anemometer 
surveys of the inlet to calculate volume flow rate. Reference velocity 
for the hot-wire reference setting is obtained from the inlet survey 
Pitot-static probe mounted at the inlet centerline. Inlet guide vane 
velocity is obtained from the permanently mounted Pitot-static probe 
ahead of the inlet guide vanes. 

This program assumes that the hot-wire instrument is set to 90 
per cent of full scale for the reference velocity at the inlet center- 
line and that the hot-wire instrument is linear. 



Input 
Card No. 


Format 


Fortran 


Description 


1 


12 


KKK 


Number of runs to be 


2 


313,2110 


M 


considered. 

Number of radial data 






MM 


points - symmetric on 
either side of centerline 
(maximum 99) 

Number of readings of total 



pressure and velocity head 
taken at the inlet centerline 
(maximum 99) 

MMM Number of readings of total 

pressure and velocity head 
taken ahead of the inlet guide 
vanes (maximum 99) 



99 



1RUN 


Run number 


ISPEED 


Rotor speed (rpm) 


3 4F10.4 TI 


Temperature in the inlet (°F) 


TW 


Temperature in test cell (°F) 


BAR 


Barometric pressure (inches 
Hg.) 


TBAR 


Temperature of mercury column 
of barometer (°F) 


4 8F10.4 PT1R 


Total pressure at the inlet 


Q1R 


centerline (inches H 0) 

2 

Velocity head at the inlet 


5 8F10.4 PT2R 


centerline (inches H^O) 
Total pressure ahead of the 
inlet guide vanes (inches 


6 8F10.4 Q2R 


h 2 0) 

Velocity head ahead of the 
inlet guide vanes (inches 


7 8F10.4 R 


u 2 0> 

Radii at which data are taken 
from the centerline outward - 
data points must be symmetric 
on either side of centerline 
(inches) 


8 8F10.4 VRL 


Values of velocity ratio of 
hot-wire anemometer to the 




left side of centerline looking 
downstream (dimensionless) 



100 



9 



8F10.4 



VRR 



Values of velocity ratio of 



hot-wire anemometer to the 
right side of centerline 
looking downstream 
(dimensionless) 



Output 

The computer progran output prints the input data in the same format 
as they are read in. 

The calculated information is presented as follows: 



Description 

Barometric pressure 

Inlet temperature 

Volume flow rate 

Mass flow rate 

Velocity at inlet centerline 

Velocity ahead of inlet guide vanes 

Mass density at inlet centerline 



CCC 



( volume flow rate) 



Units 

(lb/ft 2 ) 

<°R) 

(ft 3 /sec) 

(lbm/sec) 

(ft/sec) 

(ft/sec) 

(lbm/ft 3 ) 

(ft 3 /sec) 



(velocity ahead of inlet guide vanes) (ft/sec) 
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TABLE A-l (continued) 
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TABI£ A-2 SAMPLE OUTPUT OF PROGRAM F1X)RAT 
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APPENDIX B 



DETAILS OF FLOW PROPERTIES COMPUTER PROGRAM 

This program uses the total pressure, total temperature, velocity 
head, velocity ratio from the hot-wire anemometer, and absolute flow 
angle data obtained from the radial surveys downstream of roto;r blade 
rows and calculates flow properties based on the NISRE solution, flow 
probe velocity and hot-wire anemometer velocity. The NISRE solution 
also requires total pressure and velocity head data from the permanent 
Pitot-static probe ahead of the inlet guide vanes in. order to calculate 
mass flow rate. ; 

Input 



Card No. 


Format 


Fortran 


Description 


1 


12 


KKK 


Number of sets of data to be 








read . 


2 


3I3,7F10.3 


IRUN 


Run number 






MD 


Number of radial data points 








(maximum 50) 






N 


Number of radial points re- 








quired from NISRE solution 








(maximum 100) 






BAR 


Barometric pressure (inches 








Hg.) 






TBAR 


Temperature of mercury column 








of barometer (°F) 






RPM 


Rotor speed (rpm) 






BKG 


Assumed gross weight flow 








blockage factor (BKG ^ 1) 
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GAM 


Specific heat ratio 
(dimensionless) 


LIGAS 


Gas constant (f t-lb) /lbm-°R) 


FARB 


Arbitrary multiplication 
factor for reference axial 
velocity in NISRE solution, 
renuired to start Newton 
Raphson iteration 
(1 < FARB :£ 1,2) 


3 6F10.3 ,313 VA1 


Reference axial velocity at 
root mean square radius for 
NISRE solution, which must 
be estimated within 50 per 
cent of true value (ft/sec) 


ACC 


Accuracy to which the NISRE 
solution is to match the 

actual flow rate. 

(her 1 1 Flow rate calculated 

^ | Actual flow rate 

0< ACC .00001) 


TW 


Temperature in the test cell 
(°F) 


XI 


Inlet air temperature ( 0 F) 


R0 


Outer radius of the annulus 
(inches) 


RI 


Inner radius of the annulus 




(inches) 
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MMM 


Number of readings of 
total pressure and velocity 
head taken ahead of the in- 
let guide vanes (maximum 50) 


NREF 


The number of the data point 
at which the hot wire anemome- 

i 

ter was set to 90 per cent, 
counting from the hub outward. 


NS 


Stage number 


4 8F10.3 PTIR 


Total pressure ahead of the 
inlet guide vanes (inches 


QIR 


(h 2 o) 

Velocity head ahead of the 
inlet guide vanes (inches 
h 2 0) 



The remaining input data are from the radial surveys downstream of a 
particular rotor. Data must be entered into the program from the hub 



outward, one card per radial station. 




5 6F10.4 PTR 


Total pressure (inches H^O) 


QR 


Velocity head (inches H^O) 


TTR 


Total temperature (milli- 
volts determined from an 
iron-constantan thermocouple 
with a 32°F reference junc- 
tion) 
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ALPHA 



Absolute flow angle, taken 



positive from the axial 
direction to the direction 
of rotation (degrees) 



Output 


VR Velocity ratio from hot-wire 

anemometer (dimensionless) 

R Radius (inches) 

V 



The computer program output gives the input data in the same format 
as read in for .the purpose of checking input data for errors. 

The major part of the output is divided into three sections, the 
NISRE, flow probe, and hot-wire calculations respectively. The calcu- 
lated data for each of the three methods are presented in dimensional and 
in dimensionless form. The method used to obtain the dimensionless 
quantities is described in Section 8. 

Dimensional Form 



Name 


Description 


RADIUS 


(inches) 


PT 


Total pressure (psfa) 


P 


Static pressure (psfa) 


I? 


Total temperature (°R) 


T 


Static temperature (°R) 


V 


Absolute velocity (ft/sec) 


VA 


Axial component of absolute velocity (ft/sec) 


VU 


Peripheral component of absolute velocity (ft/sec) 


VMACH 


Mach number of absolute velocity based on calculated 
static temperature (dimensionless) 
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w 


Relative velocity (ft/sec) 


WA 


Axial component of relative velocity (ft/sec) 


WU 


Peripheral component of relative velocity (ft/sec) 


WMACH 


Mach number of relative velocity based on calculated 
static temperature (dimensionless) 


U 


Peripheral velocity (ft/sec) 


ALPHA 


Angle of absolute velocity - positive from the 
axial direction to the direction of rotation (degrees) 


BETA 


Angle of relative velocity - positive from the 
axial direction to the direction of rotation (degrees) 



Dimensionless Form 

All velocities are made dimensionless by dividing by the .peripheral tip 
velocity (U^) . All temperatures are made dimensionless by dividing by 
the compressor inlet temperature. 

R/RO Radius ratio - Radius divided by tip radius 

DPT/RHO*UT **2 (P -P )/(% p l^ 2 ) 



DP/RHO*UT **2 


(p-P a )/<% pu T 2 ) 


PHIAV 


Numerical average of the dimensionless axial 
velocity components. 
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TABLE B-l LISTING OF PROGRAM FLOW 
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SUBROUTINE HOTWlR ( NREF , VR ,GAM t PGAS , SPE E D, AL PHA , V f VA , VU , T , VMACH , U , W 
,WA, WU,WMACHt BETA,MD,R,TT,P t PT ) 

DIMENSION VR < 1 I , ALPHA ( 1 1 » V ( 1 » , VA ( I ) , VU ( I ) , T ( 1 ), VMACH < l) t U(l),W(l) t 
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TABLE B-3 SAMPLE OF DIMENSIONAL COMPUTER OUTPUT FOR NISRE PORTION OF PROGRAM FLOW 
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TABLE B-5 SAMPLE OF DIMENSIONAL COMPUTER OUTPUT FOR FLOW PROBE PORTION OF PROGRAM FLOW 
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